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Using Standard Euclidean Monte Carlo techniques, we discuss in detail the extraction of the 
glueball masses of 4-dimensional SU(3) lattice gauge theory in the Hamiltonian limit, where the 
temporal lattice spacing is zero. By taking into account the renormalization of both the anisotropy 
and the Euclidean coupling, we calculate the string tension and masses of the scalar, axial vector 
and tensor states using standard Wilson action on increasingly anisotropic lattices, and make an 
extrapolation to the Hamiltonian limit. The results are compared with estimates from various 
other Hamiltonian and Euclidean studies. We find that more accurate determination of the glueball 
masses and the mass ratios has been achieved and the results are a significant improvement upon 
previous Hamiltonian estimates. The continuum predictions are then found by extrapolation of 
results obtained from smallest values of spatial lattice spacing. For the lightest scalar, tensor and 
axial vector states we obtain masses of mQ++ = 1654 ± 83 MeV, m 2 ++ = 2272 ± 115 MeV and 
— = 2940 ± 165 MeV, respectively. These are consistent with the estimates obtained in the 
previous studies in the Euclidean limit. The consistency is a clear evidence of universality between 
Euclidean and Hamiltonian formulations. From the accuracy of our estimates, we conclude that 
the standard Euclidean Monte Carlo method is a reliable technique for obtaining results in the 
Hamiltonian version of the theory, just as in Euclidean case. 

PACS numbers: 11.15.Ha, 12.38.Gc,11.15.Me 


I. INTRODUCTION 


In lattice QCD, Monte Carlo (MC) simulations on an 
Euclidean lattice is the most popular method and is the 
preferred technique for ah initio calculations in QCD in 
the low energy regime. The calculation of the glueball 
spectrum is of major importance as it may allow an in¬ 
teresting comparison of QCD predictions derived from 
the first principles with the experimental candidates for 
the glueball states such as the /o(1710), /2(2400) or the 
77 ( 1440 ). Whilst Euclidean lattice gauge theory (LGT) 
is expected to give the most reliable estimates for QCD 
spectroscopy, the MC esults for dueball masses have still 
been an issue under debate diliiillliillilll 
There are other areas where Euclidean MC methods have 
been less successful. Examples are QCD at finite temper¬ 
ature and density, glue thermodynamics and heavy-quark 
spectra. Despite been rather neglected in comparison 
to Euclidean LGT, the Hamiltonian framework offers an 
interesting alternative [n HI HI HE m HI HI that 
needs to be explored. Hamiltonian formulation of lattice 
QCD has an appealing aspect in reducing LGT to many- 
body problem such that techniques familiar from quan¬ 
tum many-body theory and condensed matter physics, 
such as the series expansions [^, the t-expansion 
and the plaquette expansion |22l 1^ can be used to ad¬ 
dress the problem. Another advantage of this formula¬ 
tion is that from a numerical point of view, the reduction 
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in the dimensionality of the lattice from four to three 
dimensions provides a significant reduction in computa¬ 
tional overheads. 

Owing to the success of MC methods in the Euclidean 
formulation, one might expect similar levels of success 
in Hamiltonian regime. However MC approaches to the 
Hamiltonian version of QCD have been less successful 
and lag at least a decade behind the Euclidean calcula¬ 
tions. While the studies of the SU(N) glueball spectrum 
in 2-1-1 dimensions are definitely feasible, the accurate 
studies of the SU(3) glueball spectrum in 3-1-1 dimensions 
may be some way off in Hamiltonian LGT. A number of 
quantum MC methods have been applied to Hamiltonian 
LGT in the past, with somewhat mixed results. One of 
the first attempts at such a calculation was performed 
using a Green’s function MC approach pioneered by Hey 
and Stump for U(l) and SU(^ 1^ and later 

extended by Chin et al for SU(3) |^[^|^. However 
the later investigations [sE showed that this approach re¬ 
quires the use of a “trial wave function” to guide random 
walkers in the ensemble towards the preferred regions of 
configuration space. This introduces a variational ele¬ 
ment into the procedure, in that the results may exhibit 
a systematic dependence on the trial wave function. A 
projector MC approach and related “stochastic 

truncation” method [sE using a strong coupling repre¬ 
sentation for the gauge field run into difficulties for non- 
Abelian models, in that it requires Clebsch-Gordon co¬ 
efficients for SU(3) which are not even known at high 
orders. The introduction of these coefficients cause de¬ 
structive interference between the transition amplitudes 
and a version of the “minus-sign problem” rears its head 
1^. In the lack of any clear success of these methods, 
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we are forced to pursue an alternative approach. 


In our previous studies, we used standard Euclidean 
MC techniques to extract the Hamiltonian limit for the 
3-dimensional U(l) model HE 113 SU(3) lattice 

gau ge theory in 3-1-1 dimensions on anisotropic lattices 
[ 3 ^ . The idea is to measure observables on increasingly 
anisotropic lattices and then extrapolate to the Hamil¬ 
tonian limit, corresponding to At = at/os —> 0, where 
at and are lattice spacings in temporal and spatial di¬ 
rections. Applications of this method to the above prob¬ 
lems have been extremely successful 0133 and have 
given rise to the great optimism about the possibility of 
obtaining results relevant to continuum physics from MC 
simulations of lattice version of the corresponding theory. 
The objective of this work is to demonstrate the efficiency 
of this method to determine the glueball masses in the 
Hamiltonian limit of the SU(3) LGT on anisotropic lat¬ 
tices. In numerical simulations on the anisotropic lattice 
one needs the information upon the renormalization of 
anisotropy and couplings. The importance of taking into 
account the difference of scales on anisotropic lattices in 
these studies has already been discussed by Klassen , 
through the renormalization of anisotropy. In our case, 
the renormalization of Euclidean coupling Pe is also im¬ 
portant in extracting our results. This may be of some 
relevance to other studies on anisotropic lattices. We 
find that the anisotropic Euclidean results converge to 
the Hamiltonian estimates, once the difference of scales 
m has been taken into account. Specifically we find that 
this is particularly important at finite anisotropy |42| in 
the extrapolation procedure. 


The rest of the paper is organized as follows: In Sec. 
nn we briefly discuss SU(3) gauge theory as defined on 
anisotropic lattices. We also discuss our extrapolation 
procedure to the Hamiltonian limit. The details of the 
simulations, the methods used to extract the observables 
and the analysis of the data, are described in Sec. nm In 
this section we also discuss our techniques for calculating 
static potential, string tension and glueball masses from 
Wilson loop operators. The determination of hadronic 
scale To in terms of the lattice spacing using the static 
potential is also outlined. We present and discuss our 
results in Sec. im the glueball mass estimates are pre¬ 
sented; finite-volume effects are studied; results in the 
Hamiltonian limit are extracted, extrapolations of our 
Hamiltonian estimates to the continuum limit and the 
conversion of our results into the physical units are made 
in Sec. El Conclusions are given in Sec. ED 


II. ANISOTROPIC DISCRETIZATION OF SU(3) 
THEORY 

A. Action 


The SU( 3) g auge theory on the anisotropic lattice has 
the action 

SM = f E E [1 - [1 - , 

^ X iaj X i 

(1) 

where the spatial and temporal plaquettes are given by 


Pij = 3l^6Tr 
Pii = TrReTr 


Ui{x)Uj{x + i)u} {x + j)uUx) 


Ufix)U 4 {x + i)Uj{x + 4)Ul{x) , (2) 


with X being the sites of the lattice, i,j the spatial direc¬ 
tions, and Ui(x) and C /4 the link variables in spatial and 
temporal directions respectively. The couplings^ fis and 
fir are defined by 


9s 



and the anisotropy factor ^ is defined by 

r 1 as 

At at' 


( 3 ) 

( 4 ) 


Due to quantum fluctuations, the couplings and the 
anisotropy deviate from their bare values. The ^ depen¬ 
dence of the couplings gs and gr leads to an energy sum 
rule for the glueball mass, which differs in an important 
way from that which one would expect naively. In the 
weak-coupling limit, the relation between the scales of the 
couplings in Euclidean and Hamiltonian formulations has 
been determined analytically from the effective actions 
HE HE HE HE’ using the background field technique 
of Dashen and Gross HE ? Hasenfratz and Hasenfratz H3 j 
and Karsch HE obtained a mapping between the equiv¬ 
alent couplings of the Euclidean and anisotropic actions 

— = ~2" + 

9 s 9 e 

— = — + Ct(^) -I- 0{g%). (5) 

9t 9e 


For ^ = 1 (Euclidean limit), one recovers the Euclidean 
theory where gs = 9t = 9e- In the limit ^ 00 , Eqs. 

El reduce to the Hamiltonian values as obtained in Refs. 


^ We have included different couplings in spatial and temporal 
directions in order to allow the freedom to renormalize so that 
correlation lengths are equal in both directions, even though the 
lattice spacings as and at are different. 
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It is convenient to write Eq. o in a more symmetric 
way: 

= E [1 - > 

X i<.j x,i 

( 6 ) 

where 

/?« = 6/5|: V = C/?o- (7) 

In the limit ^ ^ oo, /3^ goes to the Hamiltonian coupling 
(3h- For our calculations we use 

/?? = /3-E + 3 [cs(^) + Ci-(C] + 

1 = l + ^[cs(C-Cr(C]+0(/3^^), (8) 

where (^e = ^IQe- Therefore for every {Pe,0 
there is a corresponding pair of parameters {l3^,^o)- For 
our simulations, we determine by evaluating the fac¬ 
tors Cs(C and Ci-(^) by directly calculating them in 
terms of the integrals given in Ref. ^3- The val¬ 
ues of renormalization factor r] can be determined non- 
perturbatively by matching Wilson loops in temporal and 
spatial directions]^. In our case this is however only half 
the story, as there is also the renormalization of coupling. 


B. Hamiltonian limit 

To obtain the Hamiltonian estimates we measure phys¬ 
ical quantities on increasingly anisotropic lattices, and 
make an extrapolation to the extreme anisotropic limit 
^ > oo. In the naive extrapolation procedure, one might 

assume classical values of the couplings, i.e., P = Ps = Pt 
in Eq. m and extrapolate the physical quantities to the 
^ ^ oo limit at constant p. Such a procedure is not 
strictly correct, however, at the quantum level because 
P Ps Pt due to renormalization ^ . The correct 
procedure is to extrapolate the physical observables to 
^ > oo at constant Pe- The procedure is summarized as 

follows: 

(i) Choose several anisotropies ^ for a particular Pe- 

(ii) Calculate the corresponding values of P^ and fo, 
using Eq. ®. 

(iii) Use these couplings in Eq. to do simulations 
and calculate physical observables. 

(iv) Perform a polynomial fit to the data in the inverse 

square anisotropy^ Ar^ at constant Pe and extrapolate 
measured quantities to At ^ 0 > oo ). 

As can be seen from Table Q] the difference between 
bare and renormalized couplings reaches a maximum 
around ^ = 4 (At « 0.25), where the discrepancy reaches 


^ This is because the leading discretization errors in the temporal 
direction are expected to be of order a^. 


nearly P^ — Pe ~ 0.20. Since the most anisotropic lattices 
we use are in this vicinity significant discrepancies begin 
to creep in if the correct procedure is not used. 

III. METHOD 
A. Simulations parameters 

To evaluate string tension and glueball masses, we ran 
simulations on four sets of lattices of different spatial ex¬ 
tents. One such set is listed^ in Tabled Configurations 
are generated using the mixture of Cabibbo-Marinari 
(CM) pseudo-heat bath (where one updates SU(2) 

subgroups of the SU(N) matrices) and over-relaxation 
sweeps (which involves a large change in the link matrix 
to increase the rate at which phase-space is traversed) 
with a ratio We define a compound sweep as one 

CM update followed by five over-relaxation sweeps. Con¬ 
figurations are given a hot start and then 100 compound 
sweeps in order to equilibrate. After thermalization, con¬ 
figurations are stored every 50 compound sweeps. We 
find that for physical observables, such as static-quark 
potential and glueball masses, a suitable mixture of heat 
bath and over-relaxation sweeps decorrelates field config¬ 
urations significantly, faster than a pure heat bath. The 
successive configurations were found sufficiently indepen¬ 
dent to justify the cost of evaluating glueball operators. 

B. Static quark potential and string tension 

In order to convert the lattice observables into physi¬ 
cal units, one needs to set the scale by determining the 
lattice spacing for each Pe value, where the mass of the a 
low-lying particle is typically used. Since the experimen¬ 
tal determinations of glueball masses are not clear, we use 
another purely gluonic quantity, i.e., hadronic scale pa¬ 
rameter To, instead. This quantity can be measured very 
accurately on the lattice‘s. To determine rg, in terms of 
the lattice spacing Cs, we measure the static-quark po¬ 
tential V (r) for various spatial separations f from the ex¬ 
pectation values of the Wilson loops. The Wilson loops 
are expected to behave like 

W{r,t) ='^Ai{r)exp[-Vi{r)t], (9) 

I 

where the summation is over the excited state contribu¬ 
tions to the expected value, and I = 1 corresponds to 
the contribution from the ground state. To obtain the 


^ It should be noted that Table |T|does not show all the used values 
of parameters in this set of runs in the glueball measurements. 
A disadvantage in using vq is that its physical value must be 
deduced indirectly from the experiment and there is some ambi¬ 
guity in doing this. 
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TABLE I: Run parameters used in the glueball simulations. 


Volume 

Pe 


Pi 

^0 

8® X 12 

5.5 

1.5 

5.5943 

1.4128 

8® X 16 

5.5 

2.0 

5.6379 

1.8261 

8® X 20 

5.5 

2.5 

5.6986 

2.2403 

8® X 24 

5.5 

3.0 

5.6849 

2.6683 

8® X 32 

5.5 

4.0 

5.6881 

3.5101 

8® X 12 

5.6 

1.5 

5.6443 

1.4143 

8® X 16 

5.6 

2.0 

5.7495 

1.8325 

8® X 20 

5.6 

2.5 

5.7764 

2.2525 

8® X 24 

5.6 

3.0 

5.7850 

2.6736 

8® X 32 

5.6 

4.0 

5.7882 

3.5179 

8® X 12 

5.8 

1.5 

5.8943 

1.4117 

8® X 16 

5.8 

2.0 

5.9495 

1.8378 

8® X 20 

5.8 

2.5 

5.9764 

2.2602 

8® X 24 

5.8 

3.0 

5.9850 

2.6837 

8® X 32 

5.8 

4.0 

5.9882 

3.5326 

8® X 12 

6.0 

1.5 

6.0943 

1.4197 

8® X 16 

6.0 

2.0 

6.1495 

1.8428 

8® X 20 

6.0 

2.5 

6.1764 

2.2675 

8^ X 24 

6.0 

3.0 

6.1850 

2.6933 

8® X 32 

6.0 

4.0 

6.1882 

2.5464 

8® X 12 

6.2 

1.5 

6.2943 

1.4226 

8® X 16 

6.2 

2.0 

6.3379 

1.8442 

8® X 20 

6.2 

2.5 

6.3986 

2.2669 

8® X 24 

6.2 

3.0 

6.3849 

2.7020 

8^ X 32 

6.2 

4.0 

6.3881 

3.5593 

8® X 12 

6.4 

1.5 

6.4943 

1.4245 

8® X 16 

6.4 

2.0 

6.5379 

1.8487 

8® X 20 

6.4 

2.5 

6.5986 

2.2735 

8® X 24 

6.4 

3.0 

6.5849 

2.7105 

8® X 32 

6.4 

4.0 

6.5881 

3.7167 


optimal signal-to-noise ratio, it is necessary to minimize 
the contamination from the excited states. This is done 
by using APE smearing method which is imple¬ 

mented by the iterative replacement of the original link 
by a smeared link. Following the single-link procedure, a 
particular spatial link is replaced by 


effective potential 


V (r, t) = —In 


W (r, t) 
W{r,t + 1) 


( 11 ) 


in the plateau region, while maintaining a good signal- 
to-noise ratio. A plot of the static quark potential V (T) 
as a function of separation r at Pe = 5.8 and At = 0.333 
is shown in Fig. Q] The data in this plot were obtained 
by looking for the plateau in the effective potential. 

The string tension is then extracted by fitting with the 
form 

V{r} = Vo + Kr-—, (12) 

r 


where Vq, K and Cc are fit parameters. The string ten¬ 
sion a is found according to K = ala. Only the on-axis 
potential were used to extract the string tension. The 
errors in the fit parameters were obtained using boot¬ 
strap method. We see that the data are fitted very well, 
in this instance giving K = 0.199(4), Vq = 0.682(9) and 

Cc = 0.162(1). 


2.2 

2 

1.8 

1.6 

1.2 

1 

0.8 

0.6 

0.4 

0 1 2 3 4 5 6 

r 

FIG. 1: Static quark potential V(f) as a function of separa¬ 
tion r. The solid line is a fit Eq. iini . 

The hadronic scale tq is then determined using 



^ _ / (1.65 - Be) 
Us y cra2 


(13) 


Ui(x) 


P 


{1 - a)U^{x) + (Uj{x)Ui{x + j) 

3¥=i ^ 


X u] {x + i) + u] {x)UPx - j)UPx - j) 


X Uj{x — j + i) 


( 10 ) 


where the P denotes a projection back onto SU(3) and a 
is the smearing parameter. The parameter a is fixed to 
0.7 and six iterations of the smearing process are used. 
The optimum value of smearing is found by looking for 
a value that minimizes the statistical variation of the 


Table im shows the results for ro/os along with the results 
for the string tension. Using these results, we can express 
the measured quantities in terms of tq . 


C. Glueball masses 

Our calculations are entirely conventional and follow 
essentially the methods described in Ref. nil The 
glueball estimates are obtained from the time-like corre¬ 
lations between space-like Wilson loops ^(t), 

Cit) = {^Ht)Hto)). (14) 
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where 


trices® 


m = m-{o\m\o), ( 15 ) 

$ is a gauge invariant operator capable of creating a glue- 
ball out of vacuum. At large t, and with periodic bound¬ 
ary conditions, the correlator tends to behave as® 

C(t) = Z [exp {-m(t - to)) + exp {-m {T - {t - to)))] 
-l-excited state contributions, (16) 

where m is the effective mass of the lowest glueball state 
which can be created by the operator $(t). In the present 
study, we concentrate on scalar, tensor and axial vector 
glueballs by measuring the masses in the representations 
T 2 '^, A++ and respectively. 

To determine mass m, C{t) must be calculated for large 
t such that it is well approximated by its asymptotic 
form. In principal we can use in Eq. itTCI) any oper¬ 
ator with the desired quantum numbers. However, in 
a numerical calculation the statistical errors are finite, 
and because C(t) decreases roughly exponentially in t, it 
will at large enough t, disappear into the statistical noise. 
Thus the use of a glueball operator for which C{t) attains 
its asymptotic form as quickly as possible is essential for 
extracting glueball mass. This requires the need to use 
operators that are close to the wave-function of the state 
in question. To obtain the optimal signal-to-noise ratio 
we need large smooth operators (j)i on the lattice scale. 
An optimized operator is found by exploiting link smear¬ 
ing procedure [l^, E3 and variational techniques nil 
Such operators have very good overlap with the ground 
state. 

Using a combination of single-link and double-link pro¬ 
cedures, glueball operators in each of the channels 
studied here, are constructed from smeared links C/j’(x) 
and fuzzy super-links Uj{x). A set of basic operators (j)i 
were constructed using the linear combinations of the 
gauge-invariant path ordered products of the smeared 
links about various space-like loops. In our simulations, 
these Wilson loops were measured from six smearing 
schemes, each scheme being a sequence of single-link and 
double-link mappings, making 24 basic operators in all 
in each of the channels investigated here. The optimized 
glueball operator was then found as a linear combination 
of the basic operators (j)i 


*0 

+ (to)) 


(18) 


for R = ^ (or A++), and which correspond 

to 0"'"'^, 2++ and I'' states respectively. We expect that 




I ^-m I 

= (19) 

i 


where | (pi) are the eigenstates of the Hamiltonian. At 
large temporal separation, the lowest mass domi¬ 

nates; it belongs to a glueball state which in the con¬ 
tinuum limit will have the lowest spin contained in the 
representation R. 

The coefficients in Eq. m were determined by 
minimizing the effective mass 


m^^'>{t) = - 


t-tn 


-In 




( 20 ) 


D. Analysis details 

In the analysis phase, the glueball masses were ex¬ 
tracted using a two-step procedure. First we solved the 
generalized eigenvalue problem 

C(t)v(^) = w(t)C(to)v(^\ (21) 

derived from minimization of Eq. CT . The correlation 
matrices C are then projected to the space of eigenvectors 
corresponding to the JV maximum eigenvalues with which 
satisfy the condition 


LUi(t) > £, i = (22) 

where e is a small adjustable parameter. The reduced 
correlator^ is 


C^(t) = v'^C^(t)v. (23) 


<P(t) = J^Vi(pi(t), (17) 

i 

where the coefficients vt were determined using the vari¬ 
ational technique. We computed 24 x 24 correlation ma- 


The parameter e is chosen appropriately so that we get 
rid of the unphysical modes caused by negative and very 
small eigenvalues®. We read the effective glueball mass 


® Due to the periodic boundary conditions, the correlation function 
Eq. Ilhl . is a hyperbolic cosine function of t —to. If one assumes 
that T is much larger than t — £o which we are able to evaluate 
the mass, one can neglect the contribution from exp[—m{T — [t — 
to))]. 


® Note that the vacuum subtraction is required only in the 
channel, since it has the identical quantum numbers as that of 
the vacuum. The vacuum expectation value vanishes identically 
in all other channels. 

^ The index R for the representation is now omitted. 

® We chose e in Eq. 1221 so that the reduced matrix had dimension 
between 2 and 10. The result then does not depend on e. 
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directly from the largest eigenvalue corresponding to the 
lowest energy 


meff 


= -In 


UJo{t + 1) 
UJo{t) 


(24) 


Alternatively, one can obtain the effective mass from 


'm-eff 


= —In 


'c^{t + iy 


(25) 


by setting t — to = 1 in Eff- <1^ - as Refs. [IQ 

Fig. (digraphs an effective mass showing the results of a 
single exponential to the optimized correlation function 
for the and channels for Pe = 5.8 and 

At — 0.5. The signal is seen out to the time slice 7 and 
reaches a plateau region for 1 < t < 7 in each channel. 
The data are little noisy for t > 7. It can be seen that 
two different evaluations of Eqs. and l|23) yield very 
consistent results within statistical errors. 



t 

FIG. 2: Effective mass plot for The open and 

filled symbols correspond to Eqs. 121 and 121 respectively. 


IV. SIMULATION RESULTS 

Here we focus on the lightest states because the cor¬ 
relations mediated by heavier states decrease so rapidly 
with t that it becomes hard to know whether we have 
indeed isolated the asymptotic exponential decay. More¬ 
over, glueballs that are heavy enough will decay into 
lighter glueballs and this may require more careful analy¬ 
sis. The variationally optimized glueball correllators were 
calculated on the all lattices given in Tabled For each 
glueball channel we find a region t € [tmim tmax] in which 
the correlation function was well described by its asymp¬ 
totic form Eq. To extract the masses from our 

data sets, we examine the onset of the plateaus in all ef¬ 
fective masses. In most cases, acceptable fit results were 
obtained from the fit interval t G [2,6]. 

Fig.El shows the effective scalar mass and scalar mass 
fit as a function of t for I3e=Q-0 and 6.2 at At = 0.333. 


As a result of the link-smearing and variational tech¬ 
niques, a significant signal is found and an impressive 
plateau is observed in the effective mass plot. 



FIG. 3: Effective mass plot for Af'^ channel at /Be = 6.0 
(solid triangles) and 6.2 (solid diamond) for At = 0.333. 

Fig.El shows the effective mass plot for tensor and ax¬ 
ial vector glueballs for (3e = 6.0 and At = 0.333. The 
data at lower [3e values {^e < 5.4) turned out to be 
too noisy to allow for reliable estimates of the masses, 
with no clear plateau visible and errors much larger than 
the data. A possible reason for this faster exponential 
drop-off of the correlation functions is that at the strong 
coupling (smaller Pe-, and so larger a), the statistical er¬ 
rors overwhelm the signal at smaller values of t and no 
useful effective mass beyond t = 2 is obtained. In these 
cases one simply assumes that meff{t = 2) provides a 
good mass estimate. The effective mass plateau for the 
heavier states becomes less significant due to very large 
statistical errors at smaller t values. We have not in¬ 
cluded our estimates for tensor and axial vector glueballs 
for i3e < 5.4 in our final analysis of mass gaps and mass 
ratios. The signal improves and the plateau becomes 
more visible as we walk towards the weak-coupling re¬ 
gion. One, however, does expect large finite-size correc¬ 
tion as one moves towards the weak coupling region. We 
comment on this below. We measure both the E~^~^ and 
to check the expected degeneracy and the restora¬ 
tion of rational invariance, since these states contribute 
to = 2++ in the continuum 0(3) symmetry group. 

Our results for the glueball masses in the lattice units 
for various Pe at At = 0.50 and 0.25 are shown in Ta¬ 
bles and respectively. All errors were obtained by 
jackknife method. 

Among the lattice artifacts, the finite volume effects 
pose a major problem in the glueball mass determinations 
from coarse lattice simulations. The mass of particles 
confined in a box with periodic boundary conditions can 
differ appreciably from their infinite-volume values. It is 
being suggested that finite-size effects can also induce a 
mass splitting of the E and T polarizations. 

To measure the systematic errors in our estimates from 
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FIG. 4: Effective mass plot showing the results of a single 
exponential fit of glueball correlation functions for and 
channels for Pe = 6.0 and At — 0.333. 


TABLE II: Results for the string tension K = a^a, hadronic 
scale parameter ro in terms of lattice spacing as and the 
Coulomb coupling Cc from At = 0.50 simulations. 


Pe 

Cc 

K{= ascr) 

as/To 

5.0 

0.3611(8) 

0.81(5) 

0.79(4) 

5.2 

0.341(2) 

0.73(2) 

0.75(1) 

5.4 

0.4127(6) 

0.601(7) 

0.699(7) 

5.5 

0.329(3) 

0.431(4) 

0.571(6) 

5.6 

0.296(2) 

0.273(2) 

0.449(4) 

5.8 

0.283(2) 

0.1103(8) 

0.284(2) 

6.0 

0.265(1) 

0.0683(5) 

0.222(1) 

6.2 

0.302(2) 

0.0364(2) 

0.164(1) 

6.4 

0.298(2) 

0.0227(1) 

0.129(1) 


TABLE III: Glueball masses in the lattice units from Ar = 
0.50 simulations. 


Pe 


CLtTflrpp--}' 

^2 

atm^++ 

QitTTLrpP — 

1 

5.0 

2.12(3) 




5.2 

1.74(4) 




5.4 

1.43(2) 

2.11(2) 

2.06(3) 

2.48(4) 

5.5 

1.274(8) 

1.68(3) 

1.62(4) 

1.93(3) 

5.6 

1.133(6) 

1.51(2) 

1.48(3) 

1.72(2) 

5.8 

0.87(2) 

1.11(2) 

1.075(8) 

1.35(2) 

6.0 

0.663(6) 

0.890(8) 

0.85(4) 

1.154(8) 

6.2 

0.564(4) 

0.721(3) 

0.684(4) 

0.942(5) 

6.4 

0.433(3) 

0.647(4) 

0.614(6) 

0.730(2) 


finite volume, we calculate masses on volumes OgLi and 
asL 2 which are both significantly larger than hadronic 
length scale and measure the changes in the glueball 
masses as the volume is changed. Essentially it tells us 
that if the change in masses is small when we go from Li 


TABLE IV: Glueball masses in the lattice units from At = 
0.25 simulations. 


Pe 

atm^++ 

^2 

a£?7i^++ 

1 

5.0 

2.01(1) 




5.2 

1.64(2) 




5.4 

1.427(8) 

2.10(1) 

1.99(3) 

2.37(2) 

5.5 

1.182(5) 

1.603(9) 

1.536(7) 

1.92(1) 

5.6 

1.027(3) 

1.37(1) 

1.26(2) 

1.56(1) 

5.8 

0.81(1) 

1.08(1) 

1.02(2) 

1.31(2) 

6.0 

0.598(3) 

0.841(3) 

0.80(4) 

1.116(6) 

6.2 

0.511(1) 

0.690(3) 

0.624(4) 

0.892(3) 

6.4 

0.410(1) 

0.563(3) 

0.514(4) 

0.727(3) 


to L 2 , then we can be confident that masses calculated 
on the largest volume are identical, within statistical er¬ 
rors, to the infinite volume limit masses. This strategy to 
control finite volume corrections, however, requires large 
simulations over a number of (3 e values on a range of 
lattice volumes. In practice one chooses a couple of Pe 
values to perform analysis on large volumes. This allows 
one to establish what volumes are large enough that any 
mass shift becomes visible within statistical errors. 

In order to ascertain the volume effect on our glueball 
measurements, we performed extra simulations at Pe = 
6.0,6.2 and 6.4 over a range of anisotropies, At = 0.666— 
0.25 (since we also intend to examine the Hamiltonian 
limit at these couplings). The spatial length varies from 
L = 10 to L = 16 which corresponds to a variation from 
0.90 - 1.5 fm, a range of sizes that satisfies the naive 
condition that spatial extent should be around or greater 
than 2 to. The input parameters used in these runs are 
same as those used in the 8^ x T volume simulations. The 
results from the L = 10,12 and 16 runs for the glueball 
masses in lattice units are given in Tables IVl I V II and I V III 
We plot the masses against the spatial size in Fig. El 



L 

FIG. 5: Variation of effective glueball masses with the spatial 
size L. 

There are few observations we can make from Fig. El 
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TABLE V: The glueball masses in the lattice units for various 
spatial extensions from Pe = 6.0 and Ar = 0.333 simulations. 


L 

a* 777^++ 

dt'^rp+A 

-‘2 


-'i 

8 

0.63(1) 

0.89(2) 

0.83(2) 

1.19(2) 

10 

0.69(2) 

0.99(3) 

0.94(3) 

1.26(3) 

12 

0.73(2) 

1.05(3) 

1.00(3) 

1.31(3) 

16 

0.76(3) 

1.06(4) 

1.04(4) 

1.30(4) 


and Tables El EH and Iwl 

(i) We do indeed observe substantial finite size effects, 
as our results at Pe = 6.0, 6.2 and 6.4 were found to differ 
appreciably from those obtained on larger volumes. 

(ii) The and states show a separation at 

small volume, indicating the breaking of rotational sym¬ 
metry. This separation, however, decreases as L in¬ 
creases. This has been observed for tensor glueball mass 
for SU(2) 0. The reason for this change lies in the pos¬ 
sible contamination of the small lattice glueball config¬ 
urations by those due to torelon states (flux excitations 
encircling the toroidal lattice) 0 . Thus state seems 
to behave more like a torelon state for small spatial size. 

(iii) We infer from Fig. Othat the results from larger 
volume differ very little. We find that for Pe > 6.0 a 
spatial size of L > 12 is large enough for the lightest 
glueballs to be free of finite size effects within our statis¬ 
tical errors. 


exp(-wo£) \ 

V moL } 

where pc is related to the strength of an effective triple 
scalar glueball interaction vertex. The mass shift in 
above equation arises from the exchange of scalar glue- 
balls across the boundaries of the lattice. 

Fitting the data at larger volumes to the first two lead¬ 
ing terms of Eq. yields the glueball mass in the in¬ 
finite volume limit. We did the fit first using the 
above fit function with atTOo(oo) and pc as fitting param¬ 
eters. A fit to the scalar data at Pe = 6.0 and At = 0.25 
yields atTOo++(oo) = 0.735±0.011 and po = 184±27. The 

best-fit values for atmQ++ (oo) was then used in the fits to 
the results for other two channels. A similar fit the ten¬ 
sor data gives 0 * 7712 ++ = 1.049±0.022 and p 2 = 144±19. 
We also used Eq. for the axial vector. Overall, it 
appears that at Pe = 6.0, 0*7770,1,2(8) — 0*7770,1,2(00) is 
of the order of 10 — 15%. The results of these fits are 
summarized in Table IVIIII 


TABLE VIII: The glueball estimates in the infinite volume 
limit for Ar = 0.25 at various Pe values. 


Pe 

0*777^++( 00 ) 

at?7i„++(oo) 

^9. 

atm„+- ( 00 ) 

6.0 

0.73(1) 

1.04(2) 

1.33(2) 

6.2 

0.58(1) 

0.78(3) 

1.06(4) 

6.4 

0.483(5) 

0.637(8) 

0.779(8) 


TABLE VI: Glueball masses in the lattice units for various 
spatial extensions from Pe = 6.2 and Ar = 0.25 simulations. 


L 

0*777^++ 

dlTTlrpP-A 

^2 


-'i 

8 

0.511(1) 

0.690(3) 

0.62(2) 

0.892(3) 

10 

0.570(6) 

0.747(7) 

0.709(5) 

0.974(8) 

12 

0.598(8) 

0.783(9) 

0.768(6) 

1.05(1) 

16 

0.582(1) 

0.79(1) 

0.78(2) 

1.06(1) 


TABLE VII: Glueball masses in the lattice units for various 
spatial extensions from Pe = 6.4 and Ar = 0.333 simulations. 

L 

atm^++ 

OjtTTTj rpp- A 
^2 


1 

8 

0.418(1) 

0.597(2) 

0.567(2) 

0.716(4) 

10 

0.470(4) 

0.647(5) 

0.636(5) 

0.774(6) 

12 

0.462(3) 

0.653(3) 

0.648(6) 

0.789(5) 

16 

0.467(3) 

0.643(2) 

0.641(4) 

0.775(4) 


For sufficiently large £, the scalar and tensor glueball 
masses mo{L) and 7772 (£) deviate from their infinite vol¬ 
ume limits, mo(oo) and 7712 ( 00 ), respectively, by 


atmoa{L) 


0 * 7770 , 2 ( 00 ) 


1 - PG 


exp ( — '/5moL/2^ 


moL 


The mass estimates in L ^ 00 at various anisotropies 
are then extrapolated to the Hamiltonian limit, Ar —> 0. 
To maintain the accuracy of the final results we choose 
all anisotropic points on the larger volumes and make an 
extensive analysis. This allows us to stay close to the 
Hamiltonian limit and translate to Ar ^ 0. With the 
extrapolation procedure used here, the results are almost 
as accurate at large anisotropies as in the isotropic and 
the extrapolation can be made with confidence. 


V. FINAL RESULTS FOR STRING TENSION 
AND GLUEBALL MASSES 


A. Extrapolation to the Hamiltonian limit 

The extrapolation of the string tension to the Hamil¬ 
tonian limit, performed in the powers of Ar^, is shown 
in Fig. Unlike our previous estimates |38| where we 
observed a trend towards a strong curvature in the ex¬ 
trapolation for the intermediate values of Pe, we see a 
fairly smooth variation of the string tension with Ar^. 
In the weak-coupling region K remains rather indepen¬ 
dent of Ar^, thus enabling reliable extrapolation to the 
Hamiltonian limit. 

To check the accuracy of our method, we compare our 
extrapolated estimates with Hamiltonian estimates ob¬ 
tained by using an exact linked cluster expansion (ELCE) 
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FIG. 6: Extrapolation of the string tension to the Hamil¬ 
tonian limit At —> 0. The curves are the quadratic fits in 
Ar^. 


|5f|. However, for such a comparison we must take into 
account the difference of scales between Euclidean and 
Hamiltonian regimes. The Hamiltonian coupling param¬ 
eter A = Qlgjj, where gn = lim^^oo g^ may be related to 
/3e through the relation m 

/?£; = V^- 0.07848. (27) 


We see that our results agree well with those obtained 
from ELCE although less accurate in the small Pe re¬ 
gion. We do not expect precise agreement in this region, 
since ELEC results correspond to the energy per unit 
length of a “string” of flux along one axis, whereas our 
estimates refer to the decay exponent of space-like Wil¬ 
son loops®. We do not show here the GEMC results [^ . 
which were obtained from Creutz ratios on very small 
loops, and therefore subjected to large finite-size effects. 
Our results in the Hamiltonian limit for the string ten¬ 
sion K, hadronic scale tq in terms of the lattice spacing 
Qs and the Coulomb coupling Cc are shown in Table na 
Also are shown, for comparison, the previous Hamilto¬ 
nian estimates from ELCE. 


TABLE IX: Results for the string tension K and hadronic 
scale ro and the Coulomb coupling Cc in the Hamiltonian 
limit. 


® It is well know that at the strong coupling, where rotational in¬ 
variance is broken, these two estimates d iffer at strong coupling. 
However beyond roughening transition mEi , where the ro¬ 
tational symmetry is restored, different estimates of the string 
tension should coincide. 


A 

K 

Cc 

as/ro K from ELCE 

4.2991 

0.76(3) 

0.393(4) 

0.7790 

0.664(6) 

4.6437 

0.67(2) 

0.403(4) 

0.7354 

0.563(5) 

5.0022 

0.56(1) 

0.367(3) 

0.7120 

0.404(3) 

5.1865 

0.395(9) 

0.324(3) 

0.5158 

0.34(2) 

5.3741 

0.265(4) 

0 .220(1) 

0.4311 

0.24(4) 

5.7594 

0.096(3) 

0.203(1) 

0.2580 

0.107(8) 

6.1580 

0.061(2) 

0.245(2) 

0.2090 


6.5698 

0.033(2) 

0.262(2) 

0.1542 


6.9951 

0 .021(1) 

0.251(2) 

0.1219 



To obtain estimates of glueball masses in the Hamilto¬ 
nian limit, an extrapolation of the data points obtained 
at constant Pe for various anisotropies (At = 1 — 0.25) 
is performed in powers of At®. Error estimates for the 
extrapolation may be obtained by the “linear, quadratic, 
cubic” extrapolation method |^. Figure |7| shows our 
estimates for the scalar glueballs as a function of At® for 
various fixed Pe values. Unlike the string tension, the 
scalar mass depends rather strongly on At. This can be 
seen from the slope and the curvature in the extrapola¬ 
tion to the Hamiltonian limit in the plot. This may be 
due to the presence of the 0(a®) errors in the standard 
Wilson action. In the large Pe region we, however, see a 
smooth variation of the scalar mass with At® . 



FIG. 7: Extrapolation of the scalar glueball mass (in lattice 
units) to the Hamiltonian limit At —> 0, for fixed Pe- The 
lines show quadratic fits to the data in At®. 

The extrapolation of axial vector and tensor glueball 
masses to the Hamiltonian limit is shown in Figs. Eland 
El respectively. This is performed again in powers of At®, 
in a fashion similar to the scalar glueball. In both these 
channels we note that both show a strong dependence on 
At: in the intermediate Pe region there is a factor of 2 
- 3 difference between the values at At = 0 and At = 1. 
Since the leading discretization errors in the tensor glue¬ 
ball masses are expected to be 0(a®, a®, agal), this differ¬ 
ence may be attributed to the presence of 0{aj) errors. 
However, in the extreme anisotropic region, the results 






























10 


differ very little; this is obvious since the expected 0{a^) 
errors decrease as aspect ratio ^ = 1/Ar is increased. 



FIG. 8: Extrapolation of the tensor glueball mass (in lattice 
units) to the Hamiltonian limit At —> 0, for fixed Pe- The 
lines show quadratic fits to the data in Ar^. 



FIG. 9: Extrapolation of the axial glueball mass (in lattice 
units) to the Hamiltonian limit Ar —> 0, for fixed Pe- The 
lines show quadratic fits to the data in Ar^. 


Our Hamiltonian estimates for the scalar Ms{= 
0(7710++), axial vector Ma{= atmi+^) and tensor Mt{= 
atm 2 ++) glueballs together with the Hamiltonian cou¬ 
pling A, calculated from Eq. itTTI) are shown in Table 

HI 


B. Asymptotic scaling 

An interesting question concerns the onset of “scal¬ 
ing” , which is often interpreted to mean that the physical 
quantities should become independent of the coupling in 
the weak-coupling regime. Here, we examine more closely 
the asymptotic scaling, given by the two-loop perturba- 


TABLE X: Glueball masses in the Hamiltonian limit. 


A 

Ms 

Ma 

Mt 

4.2991 

1.98(1) 



4.6437 

1.61(1) 



5.0022 

1.42(1) 

2.09(2) 

2.33(2) 

5.1865 

1.16(2) 

1.58(1) 

1.80(1) 

5.3741 

1 .01(2) 

1.40(1) 

1.79(1) 

5.7594 

0.84(1) 

1.14(2) 

1.41(2) 

6.1580 

0.738(9) 

1.034(9) 

1.32(1) 

6.5698 

0.58(1) 

0.78(1) 

1.06(1) 

6.9951 

0.47(1) 

0.64(1) 

0.79(1) 


tive /3-function 

am = C^exp^ - [l + 0{g^)] , 

(28) 

where 

. 11 . 102 

° 16^2’ ''l (16^2)2’ 

with the aim to establish that results approach the ex¬ 
pected scaling form in the weak-coupling regime. The 
coupling g refers to the Hamiltonian coupling gn in this 
case. In general, the most accurately calculated physi¬ 
cal quantity is the string tension K, and hence the di¬ 
mensionless ratio m / '/K, since the leading corrections 
to such a ratio are known to be of the order 0(1/<C)2 
0, where C is some length scale. However, it is well 
known that behavior of the string tension is quite differ¬ 
ent over the range of coupling: it follows neither two-loop 
nor three-loop scaling at these couplings and there is no 
universal Callan-Symanzik functions that can accommo¬ 
date both string tension and glueball mass down to this 
coupling range. 

Fig. ^Ishows the asymptotic scaling behavior of scalar 
glueball mass, Ms, as a function of VX. The dotted-line 
is the Hamiltonian strong-coupling expansion to the or¬ 
der /3|f and the solid line represents a fit to the Hamil¬ 
tonian weak-coupling asymptotic form It can be 

seen that our estimates appear to match nicely onto the 
strong and weak coupling expansions in their respective 
limits. The strong coupling series is not well behaved as 
the approximants begin to diverge as the crossover region 
is reached, hence allowing only rough estimates. For com¬ 
parison we also show the results obtained from plaquette 
expansion , which are found to converge quickly in the 
strong coupling region and diverge as the transition re¬ 
gion is approached. The results are somewhat lower than 
our estimates, but are almost compatible within errors. 

In the weak-coupling region, the scalar glueball mass 
appears to adhere accurately to the expected asymptotic 
scaling behavior. An unconstrained fit of the form 123) 
to the data in the range 2.3 < '/X < 2.65, which cor¬ 
responds to 5.6 < Pe < 6.4, gives a scaling slope of 
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FIG. 10: The asymptotic scaling behavior of scalar mass gap 
as a function of Hamiltonian conpling. The errors are smaller 
than the data points. 

2.83 ± 0.17 which is consistent, within errors, with the 
predicted value of 2.90. It is interesting to note that 
our estimate of Cm = 260 ± 10 agrees well with the 
value found in plaquette expansion when the curves are 
matched onto the scaling form at = 2.02. 


C. Extrapolation to the continuum limit 

To extrapolate the results given in Table IXI to zero lat¬ 
tice spacing, we combine the glueball mass estimates with 
the determinations of the hadronic scale ro/ag. Fig. 1771 
shows the dimensionless product of rp and the scalar glue- 
ball mass estimates as a function of (cs/ro)^. Also are 
shown for comparison the earlier Hamiltonian estimates 
from the series expansion ^3, t-expansion and the 
plaquette expansion |^. 

It can be seen from Fig. ^2 that scalar glueball mass 
shows a peculiar feature. There is no clear sign of the 
dimensionless quantity roTOo++ reaching to a constant 
value. As Og is increased from zero, the scalar mass first 
decreases, reaching a minimum near Ug/ro « 0.6, the 
mass then gradually increases with Ug. 

The origin of the “dip” in the scalar glueball mass is 
attributed to the sensitivity of this state to discretiza¬ 
tion errors in the standard Wilson gauge action. The 
dip is reduced by about half when classically improved 
action is used and further reduced when the classical 
improvement is supplemented by tadpole improvement 
P,@]. This implies that the origin of the dip lies, at least 
in part, in discretization errors. If this picture is cor¬ 
rect then one might expect the dip to be further reduced 
(or even eliminated ) if one accounts for the perturbative 
renormalization of action. This would remove the leading 
O(aga^) discretization errors in the tree-level, tadpole- 
improved action. However before making a definitive 
conclusion regarding the observed sensitivity, one needs 
to compute the full radiative corrections to the action. 



FIG. 11: Scalar glueball mass estimates in terms of ro as a 
function of the lattice spacing (og/ro)^. The solid and dashed 
lines represent linear and quadratic extrapolation of the data 
to the continuum limit. Results from different methods are 
also shown for comparison. 

Such calculations are under investigation and we intend 
to report on these studies in future. The fact that the 
scalar glueball is less sensitive to the finite-volume effects 
(Fig. EJ than the tensor state and somewhat larger lat¬ 
tice spacing dependence in the scalar mass could indicate 
that scalar glueball is smaller in size. 

To extract the results in the continuum limit we fit 
the mass measurements by the following three-parameter 
fitting function: 

(p(ag) = roTOG + ci(as/ro)^ -b C2(ag/ro)^, (29) 

where r^mc, ci and C2 are fit parameters. The func¬ 
tion yields a continuum limit rom.Q++ = 4.03 ± 0.05 
MeV. An extrapolation of the data, in the range 0.014 < 
(og/ro)^ < 0.18 (with 5.6 < Pe < 6.4), using a linear fit 
of the form 

(pi{as) = romo + co{as/ro)^ (30) 

yields roTOo++ = 3.94 ± 0.02. 

A comparison with previous continuum limit estimates 
obtained by Morningstar and Peardon P usiM tadpole 
improved action, Vaccarino and Weingarten [g and Te- 
per 3 using standard Wilson action on isotropic lattice 
shows that our result for the scalar mass is lower than 
that obtained by Morningstar and Peardon 0| but in 
good agreement with result obtained by Vaccarino and 
Weingarten and more accurate than the results ob¬ 
tained in Ref. 0. We, however, do not expect a com¬ 
plete agreement with the results obtained in Ref. P| 
which were computed for a different action. A compari¬ 
son with the earlier Hamiltonian results shows that our 
present estimates are hi gher than those obtained previ¬ 
ously. The t-expansion results, obtained from 1/3 
D-Pade approximants, reach a constant value near the 
continuum and are higher than the plaquette expansion 
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estimates [^, obtained by taking the average of the 
curves for various values of set parameter. A contin¬ 
uum extrapolation of the t-expansion estimates leads to 
romo++ = 3.24 ± 0.025 whereas the plaquette expansion 
estimates give 2.75±0.018. We observe that our result is 
an excellent improvement over all previous Hamiltonian 
results. 
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FIG. 12: Same as Fig. figConalpp but for the tensor glueball. 


Figs. El and El show, respectively, the continuum 
limit extrapolation of the tensor and axial vector glueball 
masses in terms of tq against lattice spacing {as/ro)^. 
Although the leading discretization errors in the ten¬ 
sor glueball mass are expected to be O{al,a 1 ,as 0 ?s), the 
tensor glueball data show less significant finite-spacing 
errors in contrast to the scalar mass. This could indi¬ 
cate that the tensor glueball has a larger size. For our 
final estimates of the tensor and axial vector glueball 
masses, we perform fits using Eq. while constrain¬ 
ing the fit parameters to be the same for both the fitting 
functions. The fits yields rom(2++) = 5.54 ± 0.06 and 
rom(l^ ) = 7.17 ± 0.17 for the tensor and axial vec¬ 
tor glueballs respectively. These results are in agreement 
with the estimates obtained in Ref. [^, however, lower 
than the results obtained in Ref. P|. This is most prob¬ 
ably due to the strong curvature in our extrapolations 
for the tensor and axial vector glueballs. In the large 
tts region our results for the tensor glueballs are in good 
agreement with plaquette expansion estimates and com¬ 
parable within errors to the t-expansion results in the 
small as region. Near the continuum our estimates are 
higher than previous Hamiltonian estimates. A linear fits 
to the t-expansion data yield continuum values 4.61±0.08 
and 5.56 ± 0.11 for roTO2++ and romi-|_respectively. 

The continuum extrapolation of the dimensionless 
mass ratio, Rt = m2++/mo++, which is expected to 
be more easily amenable, as a function of (os/ro)^ is 
shown in Figs. El The mass ratio Rt decreases slightly 
as we move from larger to smaller as and appears to 
reach a level near the continuum, but within statisti¬ 
cal errors, they can be fitted by a straight limit, giving 



FIG. 13: Same as Fig. figConalpp but for the axial vector 
glueball mass. 


Rt = 1.36 ± 0.01 at Os —> 0. In the continuum limit, our 
results agree well with the Hamiltonian results [l^l^l^ 

. This is in excellent agreement with results from the 
previous studies on Euclidean lattices. Since the Hamil¬ 
tonian formulation is thought of as a very asymmetric 
Euclidean theory, this agreement provides a test of uni¬ 
versality between the two formulations. It is interesting 
to note that mass ratio Rt in SU(2) QCD exhibits simi¬ 
lar features as that in SU(3) theory: the SU(2) Rt rises 
to a value of approximately 1.5 over the intermediate vol¬ 
ume region ^ . Our estimates for the glueball masses and 
mass ratio are shown and compared with previous results 
from various other studies in Table ED 



FIG. 14: Extrapolation of the mass ratios to the continuum 
limit. The solid line is a linear fit to the data. 


We see the evidence that the axial vector mass is def¬ 
initely greater than the scalar and tensor masses for all 
couplings analyzed here. This is in accordance with the 
claim made by Hamer M that corresponding mass ratios 
are strictly greater than one for all finite couplings. 

Finally, we convert our glueball mass estimates into 
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TABLE XL Glueball masses and the mass ratio in the con¬ 
tinuum limit. The values indicated in bold are our final con¬ 
tinuum mass estimates. The fitting function to the data is 
given in Eq. 


Reference 

Tomo++ 

rom2++ 

romi+- 

Rt 

This work 

4.03(5) 

5.54(6) 

7.2(2) 

1.36(2) 

M & P [IJ 

4.21(11) 

5.85(2) 

7.18(4) 

1.39(4) 

V & W 

4.0(1) 

5.5(2) 


1.37(5) 

Teper ^ 

3.65(11) 

5.1(3) 


1.41(7) 

UKQCD [9] 

3.78(12) 

5.53(24) 

6.90(8) 

1.44(6) 

Hamiltonian ^L| 

3.0 - 3.3 

4.25 - 4.60 

5.40 - 5.80 

1.2-1.6 


physical units by assigning value to the hadronic scale. 
The estimate for ro was obtained by combining the re¬ 
sults from various quenched lattice simulations with the 
standard Wilson action Q. The simple average = 
410 ±20 MeV of the determinations in Ref. Q was taken 
as our estimate of the hadronic scale. For the lowest- 
lying scalar glueball we obtain m(0++) = 1654 ± 83 ± 70 
MeV. The tensor and axial vector masses, for which 
the extrapolation to zero lattice spacing encountered no 
problems, are to( 2++) = 2272 ± 115 ± 125 MeV and 
to( 1'' ) = 2940 ± 165 ± 120 MeV, where the first error 
comes from uncertainty in r^mc and second error comes 
from the uncertainty in . Our estimates are in good 
agreement with those found in Ref. @ but are under the 
predicted values for the resonances /o(1710), /2(2400) by 
about 10-15%. However, for direct comparison with the 
experiment, we have to take into account the contribu¬ 
tion from corrections due to the light quark effects and 
mixing with conventional hadrons: the quenched approx¬ 
imation will easily lead to an underestimate of the sim¬ 
ulated values. A possible solution to this problem is to 
construct a QCD Hamiltonian that couples conventional 
hadrons and glueballs of identical so that physical 
states are a linear combination of hadron and glueball 
states. This will shed some light on the contents and 
dominant components of the predicted resonances in the 
scalar and excited glueball states. The glueball with ex¬ 
otic quantum numbers will, however, not mix with the 
conventional hadrons and would be ideal for establishing 
the existence of glueballs. 

VI. CONCLUSIONS 

In this work we used standard Euclidean MC meth¬ 
ods to extract the Hamiltonian limit of 4-dimensional 
SU(3) LGT. By taking the renormalization of both the 
anisotropy and couplings we have calculated the string 
tension and glueball masses in the extreme anisotropic 


limit. We have found that the renormalization of the 
couplings and anisotropy has a large influence on extrap¬ 
olating our results to the Hamiltonian limit At ^ 0. It 
has been suggested that this renormalization has a much 
smaller effect for improved actions uu , however later 
calculations showed that the discrepancy between 
the couplings /3j and I3e reach a maximum near high 
anisotropies (At < 0.333). It would be interesting to see 
how this renormalization influences the extrapolation of 
the results for improved actions to the Hamiltonian limit. 

Estimates for the scalar, tensor and axial vector glue¬ 
balls were extrapolated to the continuum limit and the 
results are presented in terms of the hadronic scale . 
Extrapolation of the scalar glueball mass to the contin¬ 
uum limit was problematic than those of tensor and axial 
vector masses. In the continuum limit the mass ratio Rt 
was observed to scale to 1.36 ± 0.2. We have demon¬ 
strated that there is a broad agreement between our re¬ 
sults and the results obtained in the Euclidean limit of 
the theory and a substantial improvement over previously 
known estimates calculated within the Hamiltonian for¬ 
mulation. This demonstrates clear evidence of the uni¬ 
versality between the Euclidean and Hamiltonian formu¬ 
lations. Our main aim for examining anisotropic lattice 
was to investigate alternative MC procedures for obtain¬ 
ing reliable results in the Hamiltonian limit, in view of 
the lack of a reliable numerical method for Hamiltonian 
LGT. We have found that standard Euclidean MC ap¬ 
proach is more successful than other quantum MC meth¬ 
ods, in particular the Greens Function MC techniques. 

From the accuracy of the results obtained here we con¬ 
clude that standard MC approach is a preferred tech¬ 
nique for calculating physical observables in the Hamil¬ 
tonian limit, just as in the Euclidean formulation. In 
order to make the Euclidean MC method a more valu¬ 
able tool in the Hamiltonian LGT, it will be crucial to 
show that it allows one to treat eventually matter fields 
in the non-Abelian case. The major disadvantage of this 
approach is the cost in computer time, as several config¬ 
uration sets must be generated to obtain one data point. 
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